Gauge-invariant Formulation of the Second-order Cosmological Perturbations 
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Gauge invariant treatments of the second order cosmological perturbation in a four dimensional 
homogeneous isotropic universe filled with the perfect fluid are completely formulated without any 
gauge fixing. We derive all components of the Einstein equations in the case where the first order 
vector and tensor modes are negligible. These equations imply that the tensor and the vector mode 
of the second order metric perturbations may be generated by the scalar-scalar mode coupling of 
the linear order perturbations as the result of the non-linear effects of the Einstein equations. 
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To clarify the relation between scenarios of the early 
universe and observational data such as the cosmic mi- 
crowave background (CMB) anisotropics, the general rel- 
ativistic cosmological linear perturbation theory has been 
developed to a high degree of sophistication during the 
last 25 years lj. Recently, the first order approximation 
of the early universe from a homogeneous isotropic one is 
revealed by the observation of the CMB by Wilkinson Mi- 
crowave Anisotropy Probe (WMAP)0 and is suggested 
that fluctuations in the early universe are adiabatic and 
Gaussian at least in the first order approximation. One 
of the next theoretical tasks is to clarify the accuracy of 
these results, for example, through the non-Gaussianity. 
To accomplish this, the second order cosmological pertur- 
bation theory is necessary. So, the perturbation theory 
beyond the linear order has been investigated by many 
authorsd S S H and is a topical subject, in particu- 
lar, to study the non-Gaussianity generated during the 
inflation!?! and that will be observed in CMB data 8]. 

In this letter, we show the gauge invariant formulation 
of the general relativistic second order cosmological per- 
turbations on Friedmann-Robertson- Walker (FRW) uni- 
verse filled with the perfect fluid. This short letter is 
prepared to show the aspect of the special importance 
of our companion paper [lCj. briefly. The formulation in 
this paper is one of the applications of the gauge invariant 
formulation of the second order perturbation theory on a 
generic background spacetime developed in Refs.|3. ITof. 
We treat all pcrturbative variables in the gauge invariant 
manner. We also derive all components of the second 
order Einstein equations of cosmological perturbations 
in terms of these gauge invariant variables without any 
gauge fixing. 



First, 



we explain the "gauge" in general relativistic 
perturbations 0, 0, Hl| . To explain this, we have to ex- 
plain what we are doing in perturbation theories, at first. 
In perturbation theories, we treat two spacetimes. One is 
the physical spacetime M which we will describe by per- 
turbations. In cosmological perturbations, M is an uni- 



verse with small inhomogeneities. Another is the back- 
ground spacetime Mo which is prepared for perturbative 
analyses. In cosmological perturbations, Mo is the FRW 
universe with the metric 



g a b = a 2 (v) {-{dvi) a (drf) b + j ij (dx ,, ) a (dx J ) b ) 



(1) 



where 7^ is the metric on a maximally symmetric 3-space 
with curvature constant K. We note that M and Mo are 
different manifolds. 

We also note that, in perturbation theories, we always 
write equations in the form 



Q("p") = Qo(p) + SQ(p). 



(2) 
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Eq. (0) gives the relation between variables on M and 
Mo- Actually, Q{"p") in the right hand side (rhs) of 
Eq. is a variable on M, while Qo(p) and 8Q(p) in 
the left hand side (lhs) of Eq. are those on A4o- 
Further, the point "p" in Q{ u p") is on M, while the 
point p in Qo(p) and SQ(p) is on Mo- Since Eq. © is 
a field equation, we implicitly identify these two points 
"p" and p through Eq. 101 by a map Mo — > M : 
p G Mo h- > "p" g M. This identification is the " gauge 
choice" in perturbation theories 

Moreover, the above gauge choice between Mo and M 
is not unique in theories with general covariance. Rather, 
Eq. @ involves the degree of freedom corresponding to 
the choice of the map X : Mo >— > M. This degree of 
freedom is the " gauge degree of freedom" in general rel- 
ativistic perturbation theory. By virture of the general 
covariance, there is no preferred coordinate system and 
we have no guiding principle to choose the identification 
map X. So, we may always change the gauge choice X. 

To develop this understanding of the "gauge", we in- 
troduce a parameter A for the perturbation and the 4+1- 
dimensional manifold Af = M x R so that Mo = Af\ x=0 
and M = M\ = A/"L_\. On this extended manifold Af, 
the gauge choice is regarded as a diffeomorphism X\ : Af 
— > Af such that X\ : Mo — » M\. Further, a gauge choice 
X\ is introduced as an exponential map with the gener- 
ator x rj a which is chosen so that its integral curve in Af 
is transverse to each M\ everywhere on Af]^. Points ly- 
ing on the same integral curve are regarded as the "same 
point" by the gauge choice X\. 
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Now, we define the first and the second order pertur- 
bation of the physical variable Q on M\ by the pulled- 
back variable X£Q on Mo, which is induced by the gauge 
choice X\ and expanded as 



X\Q\ — Qo + A £. 



£\Q 



0(A 3 ).(3) 



Qo = Q\m * s ^ ne background value of Q and all terms 
in Eq. Q are evaluated on Mq. Since Eq. © is just the 
perturbative expansion of X£Q\, we may regard that the 
first and the second order perturbations of Q are given 



by 



, and 



, respectively. 



Suppose that we have two gauge choices X\ and y\ 
with the generators x r] a and y r) a , respectively. When 
these generators have the different tangential components 
to each M\, these X\ and y\ are regarded as "'different 
gauge choices" . Further, the "gauge transformation" is 
regarded as the change of the gauge choice X\ — > ^a, 
which is given by the diffcomorphism 



A ■ = 



(4) 



The diffcomorphism $a does change the point identifica- 
tion. Further, <S>\ induces a pull-back from the represen- 
tation X Q\ in the gauge choice X\ to the representation 
y Q\ in the gauge choice 3^"a as 

y Qx = y*xQ\ Mo = (x^ytf WQ) 1 j = n x Qx- (5) 

Mo 

As discussed in Ref.Q, y Q\ — $*\ X Q\ is expanded as 
y Qx = *Q + \£s w *Q 

+y (i; 5(2) + 4 i) ) A 'Q + 0(A 3 ), (6) 

where and £? 2 % are the generators of <f>\. 

Comparing Eq. © with y% o (X^ 1 )* X Q in terms of 
the generators x r] a and y rf, we find = y rf — ^77° 
and £? 2) = [ y rj, x -n] a . Further, Eqs. © and © yield 



(2), 



( 2 ), 



q = 2£ ?(1) «g + {£ ?(2) + 4 (1) }q . (8) 



These are the gauge transformation rules for the first and 
the second order perturbations, respectively. 

Inspecting gauge transformation rules J7J) and JSJ, we 
define the gauge invariant variables of each order. First, 
we consider the metric perturbations. The metric on M.\ 
is expanded through the gauge choice X\ as 

A 2 

X\9ab = gab + X x h ab + — x l ab + 3 (X). (9) 

Since the 77 = const, hypersurfaces in the metric Q 
are maximally symmetric, the components {h^, hi V , h^} 
of the first order metric perturbation h ab is classified 



into the three sets of variables {h^, h<vL}i h(L)i ^tl}, 
{h(y) i h(TV) }> and h(TT)ij, which are defined by 



i">(VL) T n(V)i 



h(v)i> D l h( V )i — 0, 



hij =: a 2 h {L) '-f ij + a 2 h (T)ij , j^h^r)^ = 0, 
I 



l(T)ij 



= ■■ I DiDj 



A h 



(TL) 



+2D(ih(TV)j) + h{TT)iji 

D l h(TV)i = 0, D l h( TT ) i3 = 0. 



(10) 



The uniqueness of the decomposition IjlUI) is guaranteed 
by the existence of the Green functions of operators A := 
D l Di, A + 2K, and A + 3K, where Di is the covariant 
derivative associated with the metric 7^ . In terms of the 
variables given in (JTUJl, we define a vector field X a by 



X v (dr)) a + Xi(dx%, 



X n :— fyvx) - -a 2 <9 T /i( Ti ), 



(11) 
(12) 



X; := a 2 I h(rv)i + ^Dih( TL )\ , (13) 

where X a is transformed as yX a — xX a = C(i)a under the 
gauge transformation J7J). We also define the variables 

-2a 2 $ := h m -2{d n -n)X n , (14) 
-2a 2 $ := a 2 (h (L) - ^Ah (TL ^j + 2HX, V (15) 

2 (1) , 2i u 

a vi :— — a O n ri(TV)i 



(16) 

L (TT)i 3 , (17) 

where 7i = d n a/a. Further, we denote 



(1) h 
Xij ■= h( 



(i) 



Hab = ~2a $ (dr)) a (dri) b + 2a v t (drj)( a (dx l ) b ) 

W (-2 $ 7«+ Xi^ (^ l )a(rf^') fc , (18) 



where D l ^=^^.1=^.= x^ — 0- The tensor field "H a t, 
is gauge invariant under the gauge transformation (JJJ . In 

the cosmological perturbations 0, {$, v]/}, v}, Xij are 
called the scalar, vector, and tensor modes, respectively. 
In terms of the variables TL ab and X a , the original first 
order metric perturbation h ab is given by 



hab ='■ 'Hab + £ ' Xg a b 



(19) 



Since the scalar mode dominates in the early universe, 
(i) (i) 

we assume that ^i—Xij = in this letter. 

Now, we consider the decomposition of the second or- 
der metric perturbation l a b into the gauge invariant and 
variant parts. Through the above variables X a and h a b, 
we first consider the variable L a b defined by 



Lab '•= lab — %£xh a b + £ x9ab 



(20) 
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Under the gauge transformation rules J7J and (JSJ, the 
variable L a b is transformed as 

yLab - xL a b~£a9abi ^ '■— £" 2 ) + [£(1) ,^]°- (21) 

We note that this gauge transformation (|21|l have the 
same form as that of h ab - Then, as shown above, there 
exist a vector field Y a and a tensor field C a b such that 
the second order metric perturbation l a b is given by 

lab =■ Cab + 2£xh a b + {£y — £ x) 9ab- (22) 

The variables C a b and Y a are the gauge invariant and 
variant parts of l a b, respectively, and the components of 
Cab are given by 

(2) . „ (2) 



Cab 



-2a $ (df]) a (dr])b + 2a i/j (dri)( a (dx l ) b ) 



(2) 



(2) 



1-2 * X« ) (rfx 4 ) a (dx J ) 6 , (23) 

■ (2) (2) W ■ (2) 

where D l Vi—X[ij]=X l i = D % Xij= 0. The vector field 
Y a is transformed as yY a — %Y a = a a under the gauge 
transformations Q and ©■ 

Further, by using the above variables X a and Y a , we 
can find the gauge invariant variables for the perturba- 
tions of an arbitrary field as 

(1) Q := {1) Q-£xQo,, (24) 

(2) Q := {2) Q-2£ x [1) Q-{£y - £ x }Qo. (25) 

Through the gauge transformation rules J7J and JSJ, we 
can easily check that these variables are gauge invariant 
up to the first and the second order, respectively. 

As the matter contents, in this letter, we consider the 
perfect fluid whose energy-momentum tensor is given by 

f a b = (-e+p)u a u b +p5 a b . (26) 

We expand these fluid components e, p, and u a as Eq. 

'"4 



-0(A 3 



P = p + X P +^A 2 ( P +<3(A 3 ), 



(27) 
(28) 



u a = u a + A ( u a +^A 2 ( u a +0(A 3 ). (29) 

Following the definitions (|24l) and H25|) . we easily obtain 
the corresponding gauge invariant variables for these per- 
turbations of the fluid components: 

(1) (1) (1) (!) 

V-=P ~£xP, U a :=(u a ) -£xu a , 



£ 



(1) n 

e -£ X £, 



(2) (2) (1) , , , 

£ := e -2£ x e -{£ Y - £ 2 x }e, 

p := P -2£ x P -{£y-£ x }p, (30) 

(2) (2) (1 ) 

Ua := {u a ) -2£ X U a - {£y - £"x} U a . 

Through g ab u a Ub — g ab u a Ub = — 1 and neglecting the 
(i) 

rotational part in U a , we see that 

(!) (1) (1) 

U a = -a $ (drj) a +aD t v (dx l ) a - (31) 



We also expand the Einstein tensor as 

G b = G b + A«G a b + \\ 2{2) G b + 0(A 3 ). (32) 

Equations (|19|l and (|22|l give each order perturbation of 
the Einstein tensor which is decomposed as 

(1) G a b = w G a b in} + £xG a b , (33) 

{ 2) G b = wg a b[c] + ^g a b [n : n] 

+2£x W G b + {£ Y -£ 2 x }G a b (34) 

as expected from Eqs. and (J25J). Here, (1) £ a b [H] 
and ^ l 'Q a b [C] + ^Q a b [H,7{\ are gauge invariant parts 
of the frist and the second order perturbations of the 
Einstein tensor, respectively. On the other hand, the 
energy mometum tensor (|26f) is also expanded as 

f a » = T a b + A«T a b + ^ 2 >T a b + 0(A 3 ) (35) 

and ' 1 - > T' a f> and ( - 2 ' ) X , b are also given in the form 

{1)T j = iiy Ta b + £xTa ^ (36) 
(2^6 = VYr a b + 2£ x {1) T a b + {£y - £ 2 x }T a b (37) 

through the definitions (|30|> — of the gauge invariant 
variables of the fluid components. Here, ^T a b and ^T a b 
are gauge invariant part of the first and the second or- 
der perturbation of the energy momentum tensor, respec- 
tively. Then, the first and the second order perturbations 
of the Einstein equation are necessarily given in term of 
gauge invaraint variables: 

Wg a b [H] = 87rG (1) T a b , (38) 
{1) g a b [C] + ^g a b [H, H] = 8vrG ^T a b . (39) 

The traceless scalar part of the spatial component 

_ (i) (i) 

of Ecj. I|38|l yields = <J> , and the other components of 

Eq. I|38|) gives well-known equations 0. Neglecting the 

(i) (i) 

first order vector and tensor modes and using = $ . 

(2) 

the components of U a are given by 



w ( /coy (i) , (i) (2) \ , , 

U a = a $ - Di v D % V - $ (dr,) a 



f (2) (2)\ 

-a Di v + V z (dx l ) a , 



(40) 



. (2) 

where D l Vi— 0. 

All components of Eq. I|39|l are summarized as follows: 
As the scalar parts of Eq. (|39|l , we have 
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(2) (2) "if 1 ^ 

4nGa 2 £ = (-3Hd n + A + 3A - 3H 2 ) $ -r + - iA~ 1 D i DjT i ^ - -T k k 

-lnd v (A + iK)- 1 (a^dw^ - IrA , (4i) 



8vrGa 2 (e + p)A $ = -29,A $ -2WA $ +AA _1 D fe r fe - 3d, A (A + 3A)" 1 (\- 1 D i D j T i i - ^ k k ^j , (42) 

(2) (2) 1 

4irGa 2 V = {p% + 3Hd v - K + 2d v H + H 2 ) $ --A^DWjT^ 



+- {d 2 + 2Hd n ) (A + 3K)- 1 yA~ 1 D i DjT i 3 -T k k j , (43) 

(2) (2) Q / 1 \ 

*- $ = -(A + 3K)- 1 (A-iDWjTj - -T k k ) . (44) 



where 



2 2 

r := %nGa 2 (e+p)D t V A -3 A ¥ -D fe $ -3 ( d n $ ) - 8 $ A $ -12 (K + H 2 ) ( $ ] , (45) 



Ti := -WnGa 2 l£ + P I A v' +12H $ A $ -4 $ S„A $ -49, $ A (46) 

r. j := 167rGa 2 (e + p)A ( « D 3 ( v -4A $ £> 3 $ -8 $ A'A * 

+2 3A $ $ +4 $ A $ + I a, $ J +4 (2<9„W + K + H 2 ) I $ J + 8H $ <9„ $ 7/. (47) 

As the vector parts of Eq. I|39[) , we have 

8nGa 2 (e + p) v\ = \ (A + 2K) $ + - D i A.- 1 D k T h ) , (48) 



(i) (l) (i) (i) (i) „ (i) „ (i) (i) 



d v a 2 W = 2a 2 (A + 2^)- i {AA- 1 AAr fc '- Ar, fc }. (49) 



— (2) 

As the tensor parts of Eq. (|39|l . we have the evolution equation of Xij 



(d 2 + 2Hd v + 2K-A) % = 2r w - | 7<J -r fc * - 3 f AA - (A + 3A) 1 ( A- 1 D h D l T k l - \r k 



-4 (A + 2if r 1 D j) A- 1 D l D k T l k - D {i (A + 2A)" 1 £ fe r,)fc) , (50) 



Equations (|49|l and l|5U|l imply that the second order vector and tensor modes may be generated due to the scalar-scalar 
mode coupling of the first order perturbation. 

(2) 

Further, the equations i|41|l and l|44|) are reduced to the single equation for $ 

(2) 

(9 2 + 3W(1 + c 2 s )d n - c 2 s A + 2d v n + (1 + 3c 2 )(H 2 - A")) $ 

- f * +f (™) 2 + 2 # - + S (SO"} - ^» + 1) - 

-| (a 2 + (2 + 3c 2 ) H5 r) ) (A + 3A)- 1 (a^AAIV - \tA . (51) 



I 

Here, we have used the second order perturbation of the equation of state for the fluid components 

(2) (2) (2) dc 2 (I) 2 Be 2 (1)(1) &T (I) 2 , . 

V= c 2 s £ +tS +^ £ +2^ £ 5 +_ 5 ,(52) 
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(1) , (2) , ■ 
where S and S are the gauge invariant entropy per- 
turbation of the first and second order, respectively, we 
denoted that c 2 s :— dp/de and r := dp/dS. The equation 
(|51|l will be useful to discus the second order effect in 
the CMB anisotropy, for example, the non-Gaussianity 
generated by the non- linear effects 0. 

We gave all components of the second order perturba- 
tion of the Einstein equation in terms of gauge invariant 
variables. This is the main result of this paper. The simi- 
lar equations are also obtained in the case where the mat- 
ter content of the universe is a single scalar field^(|. We 
also note that we have derived the second order perturba- 
tion equations without using explicit solutions to the first 
order perturbations. Therefore, the equations i|41fl - H5U|) 
are applicable to any expansion law of the background 
universe. In particular, in the case where a(rj) — 1, these 
equations are those for the second post-Minkowski per- 
turbations. Besides the explicit forms of the second order 
Einstein equations, the gauge invariant expressions l|38|) 
and (|39|l of the first and the second order Einstein equa- 
tion are true if the decomposition formula Ijl9|l is true. 
Hence, the gauge invariant formulation of the second or- 
der perturbation shown here is also applicable not only to 
cosmological perturbations but also any other situations 
of the general relativistic perturbations. 

In the case of cosmological perturbations, Eqs. i|41|) - 
(|50|l make possible to discuss the non-linear effects in 
the evolution of the universe. In many works0, the 
second order effects have been investigated through the 
conserved quantities which corresponds to the Bardeen 
parameter in the linear theory. It will be interesting to 
clarify these conservation laws by using the gauge invari- 



ant formulation in this paper. 

Further, the rotational part of the fluid velocity in 
Eqs. H48fl of the vector mode is also important in the 
early universe because this part of the fluid velocity is 
related to the generation of the magnetic field in the 
early universe[12j. The generation of the tensor mode 
by Eq. I|5U|I is also interesting, since this is one of the 
generation process of gravitational waves. We have al- 
ready known that the fluctuations of the scalar mode 
exist in the early universe from the anisotropy of the 
CMB. Hence, the generation of the vector mode and ten- 
sor mode due to the second order perturbation will give 
the lower limit of these modes in the early universe. 

Thus, the formualtion shown here have very wide ap- 
plications and we leave these our future works. These 
issues are interesting not only from the theoretical point 
of veiw but also from the observational point of view. 
Due to its very wide applications, the gauge invariant 
formulation shown here will play a key role in general 
relativistic perturbations. 
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